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Abstract 

We discuss the inclusive jet production at HERA in the next-to-leading 
logarithm approximation. Theoretical uncertainties are considered in some 
details. We show the importance of the jet rapidity distribution to constrain 
the parton densities in the photon. A comparison is made with the recent 
HI data. 
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Several works have recently been devoted to the study of hard processes in 
photon-proton scattering. The renewed interest in such reactions is due to the 
opening up of a new kinematical domain at HERA and the publication, for 
the first time, of data on large transverse momentum jets in photoproduction 
reactions at very high energies [0-10. The importance of this process arises 
from the possibility of probing the hadronic structure of the photon in a way 
complementary to the usual studies of the deep inelastic photon structure 
function. In the latter case, one can reach the very small x domain of the 
parton densities in the photon and, in particular, directly measure the quark 
distributions whereas the gluon distribution is constrained by the evolution 
equations. On the contrary, in photoproduction reactions one directly probes 
the gluon density in the photon assuming the parton distributions in the 
proton are known from elsewhere. A quantitative study of both processes 
together should then lead to a precise determination of the photon structure 
very much in the same way as it does in purely hadronic reactions. As we 
shall see, the non perturbative input to the photon structure function will 
also be constrained by hard processes in photoproduction. 

In the following we consider the next-to-leading logarithm predictions of 
large px jet production with emphasis on the interplay between processes 
where the photon couples directly and those where it interacts through its 
parton contents. Some theoretical uncertainties are analysed. We then dis- 
cuss in which kinematical domain the future HERA data will be able to 
constrain the parton distributions. A comparison with present data from HI 
is performed and we discuss in particular the role of the gluon distribution. 
We stress the importance of the rapidity distribution in the determination of 
the photon stucture functions. 

Several theoretical papers have already appeared on the subject [[|, |J. 
Contrary to previous works we stress the phenomeno logical relevance of jet 
photoproduction data and indicate what precision the data should reach and 
what rapidity domain should be explored to extract useful physics informa- 
tion. 
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1 The photon-proton cross section 



We first consider the scattering of a real photon off a proton at a fixed center 
of mass energy y/s. The cross section for the production of a jet of transverse 
momentum px and pseudo-rapidity 77 can be written 

da <!rr L> <I<t >! ' 



dprdi] dpxdr] dpTdr] 

where the first term on the right-hand side involves the "direct" coupling of 
the photon to the hard partonic process while the second term describes the 
interaction of the photon through its " structure function" . In the next-to- 
leading logarithmic approximation these terms take the form 



da 



D 



d PTdV l=q , g 



( R ) = E / dxtFi/^M) 

-~ g 

^ ( d ° l ^ 3et + ^ KP(R- M a)) (2) 

and 

da SF r 

-( R ) = E J dx 1 dx 2 F i/p (x 1 ,M)F ]h (x 2 ,M) 

+ ^fLK?f{R-M^)\ (3) 



dprdrj 



m=<?,9 



' oM \ ( da ij -* jet a s {jj) SF/ 
2tt I \ dp T d V 2tt ij 



which exhibit the higher order corrections to the hard processes K^(R; M) 
and K^ F (R; M, /i) respectively ||. The latter term is the same as in purely 
hadronic reactions and has been calculated some time ago in |7[] while the 
former term has recently been computed and used in photon-photon colli- 
sions H . For consistency, the parton distributions in the proton, Fi/ p (x, M), 
and in the photon, Fj/ 7 (x, M), also have to be defined beyond the leading- 
logarithm approximation. The parameter R specifies the jet cone size, while 
/i and M are the renormalization and factorization scales respectively. The 
factorization decomposition eqs. involving external photons, and has 



been derived in the MS scheme || using the technics of and [TT| . 

Since one of the aims of this work is the determination of the parton 
distributions in the photon we now discuss the functions Fj/ 7 (x, M) in some 
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details. They satisfy evolution equations of type (we assume one flavor for 
simplicity) 



dFi 

with 

P®F ih = £dz pQ F ih {z) (5) 

Both the Altarelli-Parisi kernels Pjj and the inhomogeneous terms fcj have a 
perturbative expansion in terms of the strong coupling constant, 



a s (M) (Q) q 8 (M) 



*. = f^ + ^tf') (6) 



The expressions for the various terms, in the MS scheme, can be found in 
[|12]j and in flBfl , [fHH respectively. The general solution of eqs. ((D is written 
superposition of two terms, 



F ih (x, M) = F^(x, M) + F$f(z, M) (7) 

where the first term on the right hand side satisfies the full inhomogeneous 
equation and vanishes at some initial M = M value whereas the "non- 
perturbative" term obeys the homogeneous evolution equation. It will be 
discussed at length below. 

The parton distributions are defined with respect to some physical pro- 
cess which is usually taken as the deep-inelastic photon structure function 
P7(a;, Q). The relation between the structure function and the parton distri- 
butions is given by: 

F?(x,Q) = F q/j (x,Q)(l + ^M C j ] 1 \x)) + F gh (x,Q)^-C^(x) + C,(x) (8) 

The Wilson coefficients C^\x) have been calculated long time ago [15|. We 
shall always assume in the following that they are defined in the MS scheme. 
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The direct term C 7 has also been known for some time and in the MS scheme 
it takes the form 

a 4 ft 2 , /i - 1 — X 



CJx) = 6 — et (x 2 + (1 - xy) In + 8x(l - x) - 1 (9) 

Z7T V X / 

When working beyond the leading logarithm approximation we are free 
to choose the expression of C 7 (x) by modifying the factorization scheme. 
With the MS choice, as in the above equation, we observe that the ln(l —x) 
factor becomes very large near the boundary of phase space and then it does 
not appear as a correction in eq. (§) since it may become numerically larger 
that the leading terms Fj/ 7 which are enhanced by a factor l/a s . In order to 



keep the concept of a perturbative expansion useful, it is proposed in |14] to 



introduce the DIS^ factorization scheme where the choice 

C^(x)\ DISj = (10) 

is made. The parton distributions thus defined satisfy an evolution equation 
of type (fD with the inhomogeneous terms k^ replaced by 



2^«i \DIS-y - ~ F iq ® G 7 I 11 ) 



as can be immediately derived by expressing the MS parton distributions 
Fj/ 7 in terms of their DIS y counterparts 

F gh = F gh \ DISl (12) 
in eq. (HI). Note that the homogeneous terms are not affected by this trans- 



formation. It is shown in JTj] that, in this convention, the leading logarithmic 
and the beyond leading logarithmic parton distributions remain very close 
to each other over the all range in x. 

We adopt here a different approach which also absorbs the troublesome 
"large" ln(l — x) terms with the added advantages that the parton distri- 
butions we define are universal (i.e. independent of the reference process) 



and obey the MS evolution equations [16|, as all hadron structure func- 
tions in practical use today. Furthermore they are physically motivated by 
a careful analysis on how to implement the non-perturbative input. One 
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criticism which may be raised against the procedure in |TJj is the follow- 



ing. At Q = Qo the structure function F^x, Qo) is entirely given by the 
non-perturbative quark and gluon distributions which are identified with the 
quark and gluon distributions in a vector meson according to the VDM hy- 
pothesis. But there is no compelling reason to make this identification in the 
DIS 1 scheme rather than in the MS scheme (the difference is the term C 7 ) 
or even to choose another physical process where to perform the identifica- 
tion with the VDM input. A careful analysis shows that at Q = Q 
the relationship between the photon structure function and the VDM quark 
distribution should be in the MS scheme (neglecting higher order terms) 



Fj(x,Q ) = CJx) + F a N ,?(x,Q ) 



C,(x)-C (x)+F^ M (x,Q ) (13) 



with 



C (x) = 6 ((x 2 + (1 - a;) 2 ln(l - x) + 2x(l - x)\ (14) 



2tt 



In words, it means that the VDM hypothesis does not apply to F^, intro- 
duced in eq. (0), but rather to the combination F^P M = F^f + C which 
is independent of the regularization scheme. Thus, the dominant part of the 
direct term has to be absorbed in the non-perturbative input. The function 
Cq(x) arises from an analysis in the collinear approximation and it reflects 
the singularity structure associated to the external photon leg: it is therefore 
independent of the hard process in which the photon is involved. In that 
sense the definition of our VDM input is universal, independent of the con- 
sidered observable. Away from Q = Q , the function Ffj^ evolves according 
to the homogeneous Altarelli-Parisi equations in the MS scheme. 

Roughly speaking, it can be said that in the approach of [|l4j the direct 
term is entirely absorbed in the anomalous photon component with an ap- 
propriate modification of the evolution equation (DIS y scheme) whereas, in 
our approach [TB |, part of the direct term is absorbed in the non-perturbative 



input at the scale Mq and then evolves with it in the usual MS scheme. 



Throughout this work we use the parton distributions defined in |fL6 



which are adjusted to reproduce the photon structure function [T7|]-[|(J]. They 



are also in good agreement with the single jet production cross section in 
7 7 collisions at TRISTAN energies [^1|], ||. The scale M at which the 
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perturbative input vanishes is Mq ~ m 2 ~ .5 GeV 2 . The "standard set" 
is obtained when the full VDM input is used in eq. ( |13D but we have the 
flexibility of changing the overall normalization of this component so that 
its effect on physical quantities can be easily analysed. The VDM input is 
related to the pion structure functions of ||22| . 

Until now no experimental results have been presented with a fixed photon- 
proton energy All data involve a convolution on the photon energy. We 
approximate here the HI kinematical configuration |l| and convolute the 7 p 
cross section with the Weizsacker- Williams spectrum to construct the jet 
cross section in e p collisions. Taking into account the experimental tagging 
condition we use 

F 7/e (z) = -?L (l + (1 - z) 2 ) In (1 7 )g 2 L * (15) 

where m e is the electron mass and the maximum virtuality of the quasi-real 
photon is given by Q max = .01 GeV 2 . The photon momentum is further 
restricted by the condition .25 < z < .7. Let us just mention that given the 
smallness of the virtuality of the incoming photon it is a good approximation 
to neglect it altogether [^3|, [|, |24 . 



For the following studies we use for the proton structure function the 
quark and gluon distributions of ref. [p5 ]. We thus have a consistent set of 



parton distributions since the pion distributions, used in the definition of the 
VDM component, have been derived from an analysis of 7r p data using the 



parametrization of [25] for the proton. 



2 Phenomenological studies 

As is well known, all perturbatively calculated cross section suffer from scale 
ambiguities. We illustrate this point by plotting for Hera energy, y/s^ = 
295 GeV, the variation of the cross section, at fixed jet px (pr — 8 GeV/c) 
and pseudo-rapidity (rj = 0). In fig. la) the leading logarithmic predictions, 
i.e. arbitrarily setting Kf = Kf- F = 0, are shown: the plot reflects the 
monotonous variation of the coupling a s (fi) when the factorization scale M 
is fixed and the monotonous variation of the structure function with M when 
fj, is fixed. There is no prefered point on this plot. The picture in the next- 
to- leading approximation is quite different (fig. lb)) and the 2-dimensional 
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surface now exhibits a saddle point referred to as the "optimal point". This 
is seen, more quantitatively, in fig. lc) which shows the "equipotential" 
lines of the surface: in particular one observes a region of stability extending 
roughly between 2.2 < /i (GeV/c) < 4.5 and 2. < M (GeV/c) < 12. The 
same pattern holds true at other values of transverse momentum and pseudo- 
rapidity relevant for Hera energies and we find that both M and fi slightly 
increase with the transverse momentum and the rapidity although the precise 
value of the scale is not important in the neighborhood of the optimal point. 
The features of scale compensation in perturbative calculations have been 



extensively studied since the pioneering works of Stevenson and Politzer [^8 
In the case of reactions involving real photons the inhomogeneous term in eq. 
(|J) introduces an added complication [EB| which has led to some confusion in 
recent studies [0J, [27[. When calculating the higher order corrections to the 



direct sub-processes (7 q — > g q, 7 g — > q q) and requiring one of the final 
partons to be at large px one has to integrate over a kinematical configuration 
where the photon decays, for exemple, into a real anti-quark and a virtual 
quark which can be on-shell. This leads to a divergence at the pole of the 
quark propagator (collinear configuration). The factorization scale separates 
the region of high quark virtuality which contributes to the higher order 
term in eq. (0) and the region of small virtuality which is associated 

to the photon structure function and therefore to 4^-. Before QCD is 
turned on (a s = 0) in the evolution equations (U), the cross section eq. ([!]) 
is exactly independent of the scale M which is just an arbitrary parameter 
which separates the hard regime from the soft (structure function) one. When 
QCD is turned on, the InM dependence is resummed to all orders in the 
structure function while the balancing term in Kf is not, thereby leading to 

a residual M dependence in the combination + as is expected 

when working with a truncated perturbative series: the scale dependence 
is partially compensated between the higher order term of eq. (Q) and the 
leading one in eq. (|) (variation associated to the inhomogeneous term, k i} 
in eq. (01)) on the one hand, and between the two terms on the right hand 
side of eq. (§) (homogeneous term in eq. (0j)). In no case, it is expected that 
■0-r and *?5 f be separatly stable under changes of M and therefore they 

dpTdr] dpTdr] sr j o j 

cannot be physical quantities to be compared with experiment. An apparent 
stability may however be obtained for ~£ Tdrj when choosing M = /i, so that 
the variation of the coupling partially cancels that of the photon structure 
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function. 

To further test the practical effect of optimization we show in fig. 2) the 
rapidity distribution at fixed pt = 8 GeV/ c under two hypotheses: the solid 
line shows the optimal result while the dashed line shows the predictions 
based on the commonly used choice M = fi = pt- The variation does not 
exceeds 8% with a moderate rapidity dependence: the optimization decreases 
slightly the predictions in the backward region (i.e. the photon fragmentation 
region) and incresases them in the forward one. In the following, we keep 
the traditional choice M = \i = px which turns out to be quite sufficient for 
comparison with present data. Also shown in the figure is the contribution of 
the direct term, eq. (§): it is negligible except in the backward region where 
it reaches 20% of the full cross section. Needless to recall that this result is 
only qualitative as the magnitude of the direct term varies much more than 
the full cross section under changes of scales. 

We show in fig. 3a) the results of our calculation under the HI exper- 
imental conditions. The rapidity distribution is obtained integrating over 
Pt > 7 Gev/c and compared with the data . The rising rapidity spec- 
trum observed experimentally cannot be reproduced. Part of the problem 
at large r\ may be due to the contamination of the jet with particles from 
the "underlying" event which may play an important role given the rather 
low jet pt values involved. In order to test the sensitivity of our predic- 
tions to structure functions we play the following extreme games: setting the 
gluon density in the photon to gives the dashed curve with a totally wrong 
rapidity dependence while setting the gluon in the proton to essentially 
decreases the cross section without changing the shape of the curve in the 
experimentally accessible domain. We also show the predictions assuming 
a vanishing VDM input: as expected only the forward rapidity domain is 
strongly affected. The only virtue of these games is to indicate which regions 
are important to constrain the various structure functions and what precision 
the data should reach to give quantitative constraints. 

In the same spirit we plot in fig. 3b) the average momentum of the 
partons in the proton, < x p >, and in the electron, < x e >, as functions of 
the rapidity, for the cross section shown in fig. 3a). We note the small value 
of < x p > over the whole range confirming the important role played by 
the gluon in the proton as just seen above. When Hera data become precise 
enough they can effectively be used to constrain the proton structure function 
in the vicinity of x p ~ .05 at values of Q 2 around 50 to 100 GeV 2 . The average 
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value < x e > of the partons in the electron is obtained by convoluting the 



photon spectrum eq. (jT5j) with the parton distribution in the photon and it 
is calculated here using only the SF component. The rather "large" < x e > 
obtained is due to the fact that the parton distributions in the photon are 
not as peaked as that in the proton. In the range considered by HI the cross 
section is mostly sensitive to the quark density which is also probed in studies 
of F2 . To achieve good sensitivity to the gluon distribution would require a 
larger rapidity coverage in the forward region. 

In fig. 4), we plot the jet distribution in after integration over the 
rapidity domain —1. < rj < 1.5. Excellent agreement is reached both in 
magnitude and in shape. In view of fig. 3a) however it is clear that this 
agreement is purely fortuitous! This illustrates the fact that rapidity distri- 
butions are more powerful than transverse momentum ones to constrain the 
parameters of the theory. 

In conclusion, we find the study of jet photoproduction at Hera very 
promising. Concerning the quantitative comparison between theory and ex- 
periment it appears that the relatively low values of the jet transverse mo- 
mentum do not allow us to exclude the possibility of contamination of the 
jets with low pt fragments of the proton making the extraction of the par- 
ton momentum from the jet momentum very difficult and the comparison 
with the theory rather dubious at this point of the analysis. This shows also 
the relevance of collecting data on inclusive particle production || at large 
transverse momentum which do not suffer from such an ambiguity. The ex- 
tra degrees of freedom contained in the fragmentation functions are by now 
rather well constrained and the comparaison with theoretical calcula- 
tions |3C] on particle production would also be very useful to determine the 
parton distributions in the photon. 
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Figure Captions 



Fig. 1 Variation of the cross section at fixed pr = 8 GeV j c and rj — 

as a function of the factorization scale M and the renormalization scale /i: 

a) leading logarithmic approximation; b) next-to-leading logarithmic 
approximation; c) contour plot for the NTL approximation. The equal 
cross section lines are labeled in pb (GeV / c)~ r ) . 

Fig. 2 Rapidity dependence of the cross section at fixed pt = 8 GeV/ c: the solid 
curve is the optimized prediction while the dashed curve is obtained with 
the choice M = \i = pr- The dotted curve is the contribution of the direct 
term (multiplied by a factor 4). 

Fig. 3 a) Rapidity dependence of the theoretical predictions compared to the 
experimental data of HI. The dashed curve is obtained when F 5 / 7 = 
and the dotted one when F g / p = 0. The dash-dotted curve results from 
the choice F^ M = in eq. ([13]). 

b) Average values of the scaled longitudinal momenta in the proton and 
the electron. 

Fig. 4 Transverse momentum dependence of the cross section integrated over 
the pseudo-rapidity range —1. < 77 < 1.5 and comparison with HI data. 
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